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В работе рассматривается осесимметричная контактная задача 
теории упругости о вдавливании жесткого сферического 
штампа с неровной формой поверхности в трансверсально 
изотропное полупространство с функционально градиентным 
трансверсально изотропным покрытием. Модули упругости в 
покрытии изменяются с глубиной по произвольным непре- 
рывным или кусочно-постоянным независимым друг от друга 
законам. Решение задачи сведено к решению парного инте- 
грального уравнения, шероховатая структура штампа модели- 
руется отрезком ряда Фурье-Бесселя. Трансформанта ядра 
аппроксимируется произведением  дробно-квадратичных 
функций, для которых получены замкнутые аналитические 
решения парных интегральных уравнений. Полученные ре- 
шения асимптотически точны для малых и больших значений 
относительной толщины покрытия. Разработана схема чис- 
ленного построения функций податливости для случая одно- 
временного действия произвольных осесимметричных нор- 
мальной и касательной нагрузок. 
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шероховатость, покрытие, неоднородность, функционально 
градиентные и слоистые материалы. 
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Введение. Контактным задачам для функционально-градиентных материалов и покрытий посвящено большое 


количество современных исследований. Однако, большинство известных в литературе результатов получены лишь 


для частных случаев изменения упругих свойств в покрытии. Возникающие при этом интегральные уравнения реша- 


ются в основном численно и полученные решения эффективны лишь в некотором ограниченном диапазоне значений 
геометрического параметра задачи (относительная толщина покрытия) [1-3]. 


При проведении эксперимента по наноиндентированию используют штампы с различной формой наконечни- 


ка. Наиболее распространены штампы сферической, пирамидальной или конической формы. Однако, в процессе изго- 


товления штампа не удаётся добиться идеальной формы наконечника, кроме того форма штампа может изменяться в 


процессе эксплуатации [4]. Перед проведением экспериментов форма наконечника индентора обычно отдельно изуча- 
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ется с помощью, например, растрового электронного микроскопа [5, 6]. Однако, математическое моделирование, поз- 
воляющее учесть отличие формы штампа от идеальной, изучено очень слабо. 

В настоящей работе разработана модель, позволяющая учесть неровности поверхности штампа и непрерывно- 
неоднородную или кусочно-однородную анизотропную структуру покрытия. Рассмотрена контактная задача о вдав- 
ливании недеформируемого штампа в трансверсально-изотропное полупространство с функционально-градиентным 
покрытием. Основное отличие от аналогичной задачи для изотропных материалов состоит в схеме построения транс- 
форманты ядра интегрального уравнения (функции податливости среды). В работе предложена схема построения 
функций податливости при действии произвольных осесимметричной нормальной и касательной нагрузок. Считается, 
что форма штампа неидеальна, неровности поверхности модулируются отрезком ряда Фурье-Бесселя. Построено при- 
ближенное аналитическое решение контактной задачи, асимптотически точное для малых и больших значений отно- 
сительной толщины покрытия и обладающее высокой точностью для покрытий средней толщины. 


Построение функций податливости. Рассмотрим упругое неоднородное полупространство © с верхней гра- 
нью Г. С полупространством связана цилиндрическая система координат г, ф, <. Ось < нормальна поверхности Г и сов- 
падает с осью анизотропии. Модули упругости полупространства изменяются по законам: 


& ©) -Н<:<0 


т — ‚ (КЛ) = 11 12, 13, 33, 44, 
и сб? = сопя —©<5<-Н ь 


где с (э) — непрерывно-дифференцируемые функции, определяющие закон изменения упругих модулей в покры- 


тии (—Н <<<0), ен" (2) — постоянные, определяющие значения упругих модулей подложки. Здесь и далее индексы 


(с) И (5) соответствуют покрытию и подложке соответственно. Покрытие и подложка жестко сцеплены между собой: 
(с) (5) (5) — <) © — 2) 
‚0: 20, ‚т =Ты. (р 


Рассмотрим действие произвольной осесимметричной нормальной и касательной нагрузок в круговой области 
0<7<а покрытия. Вне этой области поверхность не нагружена: 


&=-Н: и® =и®, и® =и 











_ |- Р@),г<а _ | (г), г<а 
25=0 | 0, г>а ’ "== |0, га’ (2 
Определяющие соотношения для трансверсально-изотропного материала имеют вид: 
с, еее", в о еыЫ, 
дг т = ? дг г 0 ) 
а с] о. ны ОИ. и с] Е "| 
" де г ` 0 . 0% ог 
Уравнения равновесия ввиду симметричности относительно координаты ф принимают вид: 
о ЗАВ =, Е Е 0. (4) 
ди 05 т д 05 т 
Используя преобразование Ханкеля 
и(г, =) = Е. и(у, <) (уау, Ди, <), (г), т} = | {Су 5), РОС} осу (5) 


и соотношения (3), уравнения (4) можно преобразовать в систему двух обыкновенных дифференциальных уравнений 
второго порядка с переменными коэффициентами: 





— , п и! 2 _ 
У (сз + су + еду +" си +и' си фисиу” =0, 








Г 7’ 2 А РИ (о 
и"сзз Ни’сзз — ИСУ” — и (св + —исз1 =0. 
Граничные условия (1) и (2), используя (5), принимают вид: 
2=-Н: © = ® п® = п, (1) 
&=-Н:сФу® ся = св Фи, (8) 
ИЕ: с@ © 4. у<©)) — с. (и + у‘) ), (9) 
5=0: со оби =-р(7), ё® и +9) = *(у) (10) 


Перепишем систему (6) в матричном виде: 


х = АС). ,-Н<5<0, (11) 
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0 1 0 0 
(с) (с) (с) (с) (с) т 
2С 44 С44 3 +544 и 
(с) (с) (с) (с) т" 
А(с) — Сад С4д С4д С44 Ее и | 
0 0 0 1 в м 
(с) (с) (с) (с) (с) 
СЗ Ч 44 2 С44 33 тя 
(с) (с) (с) (с) 
Сз3 Сзз Сзз Сзз 


Для однородной подложки (< <-Н ) система (11) упрощается до системы обыкновенных дифференциальных уравне- 
ний с постоянными коэффициентами: 








= А). х5<-Н, (12) 
0 1 0 0 
(5) (5) (5) 
с сз’ +С 
у? 11 0 0 — 13 44 
с($) с(5) 
А($) = 44 44 
0 0 0 1 
о рае в | 
(5) (5) 
Сзз Сзз 


Помимо граничных условий (7)-(10) полагаем, что выполнены условия затухания смещений на бесконечности: 
ш(г,<) —> 0, их, <) > 0. (13) 
3—— => 
Будем искать решение систем (11) и (12) в виде: 
Хз) = а“ (у, з + (уу аз (у, =) 4 
Жо (уз)= РУ (уз) (у у а* (у, =). 


Тогда, ввиду линейности, получаем системы для определения векторов а в Е 


(с) (с). (С) + 
а" = А‘ а’, 1=12, (15) 
(5) — 4(5).1(5) + 
а’ * =.“ а’, 1=12. (16) 
Общее решение системы (16), соответствующей подложке, имеет вид: 
(5) ы 9.5 
а; (1,2) = У Ри(уие . (17) 
Здесь К ж ЕК, ( 1=12; К = 1,2,3,4) — известные постоянные, зависящие от упругих модулей подложки, их значения не 


приводятся ввиду громоздкости; 0 — корни биквадратного характеристического уравнения: 











2 2 2 
а (сисзз— сз — 2азс4д )+ \ С13 — СИ6з3 (4сла (сз + сад) + с — сисзз ) 
к== . 
233С4д 


> 2 
В общем случае а, различны. Кратные корни возникают, если выполнено одно из условии: С|3 =С11Сзз ИЛИ 


2 р 
(сз +244)” = си1Сзз ‚ например, в случае изотропного материала. Подробно процесс построения функций податливо- 


сти для случая кратных корней описан в работе [7], поэтому далее будем полагать, что а, различны. Параметры О;(у) 
подлежат определению из граничных условий. Из (13) следует, что коэффициенты при слагаемых, соответствующих 
а; с отрицательной вещественной частью, равны нулю. Обозначим корни с неотрицательной вещественной частью 
через о и 9>. Тогда (17) принимает вид: 


а й 7,2)= в + Ру (уЕ ре*2К. 


Граничные условия относительно векторов а(°),а() принимают вид: 
И] 1 


<=-Н: в = =а(”, 7=1,2, (18) 


2=-Н: Фа - Фа = а Фа, (19) 


&=-Н: са ++ а) =с® м. ++ заб ), (20) 
(с) (с) (с) 1(с) 
п. | бзз ад -@вз Зал | _ _ [1 _ [0 21 
О ее, | ео | - 
Сад (ар + 7а ;з ) 
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Таким образом, для вычисления неизвестных а = [2 а’ И аа? получены две краевые двухточечные задачи 


для систем обыкновенных дифференциальных уравнений с переменными коэффициентами (15) и граничными усло- 
виями (19)—(21). Неизвестные О;(у), ],К=1,2, находятся из условия (18). Введем функции эффективных упругих моду- 
лей: 








№ С44 мСиСзз — 613 
9›(=)= ни + , 
^Сзз \Усисзз + сз +2644 


д/с А/сисзз = С 

44 1133 — С13 

®.(=)= и +ез , 
УС \Усисзз +1 +2644 


9,.(&)=©9, (<) =-усисзз - аз 


и отдельно введём обозначения для их значений, соответствующих Подложке и поверхности покрытия: 








ео =©,(0), 9 = 9, (-Н —0). Наконец, введём функции податливости аналогично терминологии А.К. Приварнико- 
ва [8]: 

[у (у,=)=@( а (у,=). 
Функции [4, (у. =), К=1,2, /=1,3 не зависят от приложенной нагрузки. В общем случае они могут быть вычислены лишь 


численно при фиксированном значении %. При ==0 функции податливости положительны для любого у, и выполнено 


Шт 1; (7,0) =1 и вм 1, (5,0) = {ый / © . Таким образом, значение функции податливости в нуле — есть отношение 
> у—0 ` ы Ч 


значений эффективных упругих модулей, соответствующих поверхности покрытия и подложке. 
Из (14) получим формулы для смещений точек покрытия в виде линейной комбинации образов преобразова- 
ния Ханкеля от действующих на поверхности нормальной и касательной нагрузок: 


ыы | Г, |= 
и (у, )= 13) 5) о 5) 














911 795’ 
=) __ 32). 1532). 
Е Ру) + во 5) (22) 


Описанная выше схема построения функций податливости может быть использована для сведения смешан- 
ных задач теории упругости к решению интегральных уравнений, при этом [„,(у,0) будут являться трансформантами 


ядер этих интегральных уравнений. 
Задача о вдавливании штампа. Пусть недеформируемый сферический штамп с неровной поверхностью 
контактирует с верхней гранью Г полупространства © по области < =0,г<а. К штампу приложена вдавливающая 


сила Р, ось которой совпадает с осью <. Силы трения предполагаются отсутствующими [рис. 1]. Под действием силы 
Р штамп переместится в направлении оси < на величину -6. Считаем, что распределение неровностей (шероховатости) 
штампа не зависит от угловой координаты и может быть смоделировано отрезком ряда Фурье-Бесселя. Требуется 


определить распределение контактных нормальных напряжений под штампом: с. в—0 = —Ра( Г ), г<а. 








Рис. 1. Постановка контактной задачи о вдавливании 
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Граничные условия при данной постановке имеют вид: 


57 =0, г>а, 


&=0: (с) = 0, з м (23) 
м2) =-б+к [26+ УВ Ло( ик г/а), г<а, 
К 


где ^ — радиус сферического штампа (который предполагается существенно больше, чем область контакта, и аппрок- 
симируется параболоидом вращения), их — положительные нули функции Бесселя Ло. 
Используя (3), (5), (22) и (23) запишем интегральное уравнение задачи: 


1 со ’ 
Грот 0] [шк 


^9(57 2..2 м 
б — Хр Лики) г’ <1 
= 








(24) 





а 


Здесь ^=Н/а, "=т/а , р(х)= ра(ха). 
Используя двухсторонний асимптотический метод [9, 10] получим приближенное аналитическое решение ин- 
тегрального уравнения (24): 


(5) 
и м И +УСАА АЯ а 


ия Г м Ил у 12 
Г г (25) 
29 м Ви, | эта) 


па Иру) и? 

Постоянные С определяются из системы линейных алгебраических уравнений: 
— З 
Ги(0)К м Ь, #9: вы (ВЕЛ! +1 


а О0<г <1. 











>в "ЕВ, А мы )+ 
7: | 








24? (ми,) СА” В} 
Здесь К=1,...,№, т — мнимая единица, Г\(и) — А трансформанты ядра интегрального уравнения: 
2; 
Ем(и)= в — (т 
ТИ В 
Вс и)+н5й(н) _ сии) 
Е(В,н)= — | 
В“-и В 


Соотношения, связывающие радиус сферического штампа, зону контакта, осадку штампа и вдавливающую силу име- 
ЮТ ВИД: 


) 


у с05( му ) 


2 (м м 
Е [$ и вт) 
К = 2 = К= РМ, ) 


в> 8а`013) м 


1+35С(сЩАХ" )- АРАЯЙ( АЙ | 
1=1 





+40(а у мь (со нуну н) )) 
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Полученные формулы являются асимптотически точными при /-—0 и ^—о0 [9]. Алгоритм построения аппрок- 
симаций высокой точности и связь между погрешностью решения для произвольного значения /. и погрешностью ап- 
проксимации трансформанты ядра интегрального уравнения описаны в работе [11]. 


Выводы. Исследования, представленные в работе, являются естественным продолжением полученных ранее 
результатов для случая изотропных материалов [12—14]. Для штампов с плоским основанием или конической формы 
функции податливости среды остаются неизменными, изменяется лишь правая часть интегрального уравнения (35). 
Решения соответствующих интегральных уравнений построены авторами ранее [7, 13]. 

Схема построения функций податливости, предложенная в работе, может быть использована при решении ря- 
да контактных задач с учётом одновременного действия на поверхности нормальных и касательных нагрузок и задач с 
учётом сил трения. Также аналогичная схема может быть применена при решении контактных задач в рамках элект- 
роупругости, термоупругости, магнитоупругости и т.д. 
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Метод решения интегрального уравнения, используемый В работе, позволяет с высокой точностью построить 
решения контактных задач даже при сложном немонотонном изменении упругих СВОЙСТВ В покрытии [15], а также в 
случае, когда упругие модули подложки более чем на порядок отличаются от упругих модулей покрытия [12]. 
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